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NUMERICAL CALCULATION OF TRANSONIC BOATTAIL FLOW 
by Wen L. Chow,* Lawrence J. Bober, and Bernhard H. Anderson 

Lewis Research Center 


SUMMARY 

A viscid-inviscid interaction procedure for the calculation of subsonic and transonic 
flow over a boattail has been developed. This method couples a finite -difference invis- 
cid analysis with an integral boundary -layer technique. Results indicate that the effect 
of the boundary layer is as important as an accurate inviscid method for this type of 
flow. Theoretical results from the solution of the full transonic -potential equation, in- 
cluding boundary layer effects, agree well with the experimental pressure distribution 
for a boattail. Use of the small disturbance transonic potential equation yielded results 
which did not agree well with the experimental results even when boundary layer effects 
were included in the calculations. The strong interaction characteristics of this type of 
problem are fully illustrated in this study. 


INTRODUCTION 

To determine the performance of an air-breathing propulsive system for a super- 
sonic aircraft operating at subsonic cruise conditions, it is necessary that the flow field 
associated with transonic flow past boattails be successfully analyzed. This type of 
problem is extremely difficult as the governing inviscid-flow equation is of the mixed 
type. Furthermore, the relatively short boattail is usually immersed within the thick 
viscous layer associated with the upstream flow so that the coupling interaction between 
the viscid and inviscid streams cannot be disregarded. Experimental simulation of such 
flow conditions has been carried out, for example, by Shrewsbury (ref. 1). Indeed, the 
experimental data indicated considerable influence of the boattail juncture shape on the 
pressure distribution on the afterbody, especially at high subsonic free-stream Mach 
numbers which, in turn, illustrated the extremely sensitive character of the transonic 
flow. 

* Professor of Mechanical and Industrial Engineering, University of Illinois, Urbana, 
Illinois. 



Since the numerical treatment of external transonic flow using small disturbance 
theory was made by Murraan and Cole (ref. 2), there has been a considerable amount of 
activity in this area. Krupp and Murman (ref. 3) computed transonic flow past lifting 
airfoils and slender bodies. Bailey (ref. 4) also extended the small disturbance treat- 
ment to calculate transonic flow past slender bodies of revolution. Steger and Lomax 
(ref. 5) and South and Jameson (ref. 6) used the full potential equation to calculate the 
transonic flow past two-dimensional and axisymmetric configurations, respectively. 
However, none of these methods considered the effects of boundary- layer growth on the 
bodies. 

The present study is restricted to the transonic flow field associated with a particu- 
lar configuration (see fig. 1) used in the experimental program reported in reference 1. 
Because of the large length to diameter ratio of this configuration, small disturbance 
theory should be sufficient to calculate the boattail pressure distribution. Such an anal- 
ysis has indeed produced accurate results for the flow past a parabolic arc of revolution 
at transonic speeds (ref. 4). These results were reproduced during the present study. 
This method was applied to the boattail configuration shown in figure 1, and the results 
are shown in figure 2. The results of the inviscid calculation do not agree well with the 
data, and the results of the small disturbance theory with boundary layer displacement 
effects, although in better agreement with the data, do not adequately predict the pressure 
distribution. This lack of agreement is probably due to the relatively large local curva- 
tures and body slopes present on the body causing large disturbances in the flow field. 

It is obvious that the small disturbance treatment is not adequate for this boattail 
configuration and the full potential equations must be used to study these problems. The 
results obtained from the present study will illustrate the strong interaction character of 
these problems within the transonic -flow regime even though flow separation has not been 
indicated on the boattail. 


FUNDAMENTAL EQUATIONS 


Inviscid Flow 


It has been well established (ref. 7) that, even with the presence of shock waves, 
the axisymmetric, inviscid, transonic flow field may be described by the potential equa- 
tion that may be written as (see appendix for definitions of symbols) 
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where the subscripts indicate partial differentiation with respect to r and z and where 
c is the local speed of sound obtained from 


c 


2 




Upon defining 


* = V co (z + <p) U = u/V w V = v/V M C = c/V^ 
the normalized disturbance potential cp satisfies 



where 


and 


u = 1 +< ?z 

V = <P r 

C 2 = — + 2UlJ: (1 - u 2 - V 2 ) 

™2 2 


Upon introducing the transformation (ref. 6) 


( 2 ) 


(3) 


(4) 


(5) 


£ = Z ? = — 2zz_ (6) 

1 + B?] 

where 

V = r - r b (z) (7) 

r b (z) is the local radius of the body, and B is a stretching parameter, equation (4) be- 
comes 
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In the £ -coordinate system the normalized velocities are given by 


and 


U = 1 + <p^ - r^B(l - C ) 2 ^ 

► 

V = <p ? B(l - S ) 2 


(9) 


where r^ and r^’ are, respectively, the first and second derivatives of r^(z). This 
sheared cylindrical coordinate system is used throughout the flow field since boattail 
models are usually provided with pointed rather than blunt forebodies. Thus, an orthog- 
onal coordinate system, as used by South and Jameson (ref. 6) for blunt forebodies to 
insure numerical stability, is not needed. The flow field to be considered in the ij, £ - 
plane is now -°° < £ < +°°, 0 < £ < 1. It is the experience derived from this study that, 
with the body situated from zero to one along the axis, the limits of £ = -2 and £ = 3 
can be treated as the upstream and downstream limits in the | direction. The physical 
and computational planes are shown in figure 3. 

Physically, the flow should be undisturbed at very large distances from the body. 
Thus, the disturbed potential should vanish at £ = -2, £ = 3, and also at £ = 1 (?j — °°). 
The boundary condition on the body surface is given by 


V - r^U at £ = 0 (t? = 0) 


( 10 ) 


which may be reduced to 


C = 0) = 


4 + ■t, 2 ) 


'l + C = 0) 


(11) 
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For locations along the £ -axis where £ = 0 and r^ = 0 (upstream of the body and 
downstream of the body when the model has no sting support) the term cp r /r in equa- 
tion (4) is indeterminate since both cp r and r are zero. In the limit as r — 0 the 
term cp r /r — cp rr and the terms in equation (8) are correspondingly modified. For 
convenience, equation (8) is rewritten as 

EF<p££ + + S = 0 (12) 

where 


E = B 2 (l - S) 4 




and 

S = </> ? B(l - 0 2 


Equation (12) can be solved by a numerical relaxation scheme. It was pointed out by 
Murman and Cole (ref. 2) that in solving the transonic flow problem by numerical calcu- 
lations, different finite-difference formulations should be adopted depending whether the 
local flow is subsonic or supersonic. For the present problem variable grid size is 
used in formulating the finite differences in the £ -direction, and uniform size is 
adopted in the £ -direction. For the grid point situated at and the equivalent 
finite difference forms of <p <p <p^, (p^, and <p ^ may be written as 


r b + 


L 

b(i - o 


- fl - — 1 - 2B(1 - OF 

c 2 
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(Pt = <pf * + — - = cpf - (0. 1 

5 ^^ 2 (A ^l +A ^ 2 ) J A ^1 A ^2 J J_1 


= ABl • ^ +1 + AB2 • <p* + AB3 • cp^ 

*« - <pf +1 <pf + — -* V* 


11 A| 2 (A4 1 + A{ 2 ) j+1 A4 jA 4 2 j A^(A| 1+ A| 2 ) 

= BB1 • + BB2 • cp* + BB3 • 
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for locally subsonic flow, 
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(15) 


for locally supersonic flow, and 
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where 


4 *2 - «j + l ' «j 


4 *1 * { i - «J-1 


A io - ?j-i - 5j-2 

Equation (12) can be written in finite-difference form for the point (ijj, and is given 
by 


E- F + 


AB2^j C+1 + (-2E- F + A? 2 D • BB2)<?jf + ^E- F - 

. - („jft - ,»$) * AB3 ■ („$ - )] 


[d.(, 


A£ 2 p- (BBl • 1 + BB3 ■ tfi 


f-i) + s ] 


for subsonic flow, and 


^E> F + APl^* 1 + (-2E- F + A? 2 D ■ BPl)^ + ^E- F - ■ APlj of ' 1 

* ■ i Y i '[ AP2 '(' p M " ’’f-i) +AM - (^2-^-2)] 

- At; 2 jp • (BP2 • + BP3 * ^. 2 ) + S ] ( 19 > 


for supersonic flow. The boundary condition (eq. (11)) is introduced for all points on 
the 4 -axis (C = 0, k = 1 for all grid points on the (j-axis) by introducing a row of grid 
points at ? = -A£ whose values of <p (identified as y ? ) are given by 

<py - ^ =2 - 2 A?<P ? Uj» K = °) (2°) 

where <P^(Ky ? = 0) is computed from equation (11) with <p ^ again evaluated depending 
whether the flow’ there is supersonic or subsonic. 
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Since the line relaxation process will be used for the numerical calculations, equa- 
tion (18) or (19) may be written for all grid points in the same column with the result 
that the coefficient matrix for the system of equations of <p is tridiagonal. These equa- 
tions may be easily solved by standard and efficient methods. The information at a grid 
point used to account for the boundary condition (£ = -A?) is updated immediately after 
the calculations for this column are completed. 


Viscous Flow 

It was noted previously that the viscous layer will be expected to modify and in- 
fluence significantly the inviscid flow field; an accurate description of the growth of the 
turbulent boundary layer is thus necessary. Although more sophisticated methods of 
prediction of its growth are available (refs. 8 and 9), a simple calculation of the growth 
of the turbulent boundary layer by the integral method of Sasman and Cresci (ref. 10) 
and then modifying the body into an equivalent inviscid geometry by adding the displace- 
ment thickness to the original body are sufficient to account for viscous effects in the 
flow. Essentially, a pair of ordinary differential equations describing the growth of 
the momentum thickness 8 and the form factor H of the boundary layer may be 

C 1/ 

written as (ref. 10) 


^ = 1.268 
dx 


t dM , dr. 

— -(l + g w Hj) - — — - + A| 

M g dx w 1 r^ dx 1 


( 21 ) 


and 


dH. 

l 

dx 


1 ^e 


2M e dx 


[h^h. + i) 2 (h. - 1)] 1 


Hf + 4H. - 1 
1 + <g„ " 1) 1 


w (H. + lXHj + 3)_ 
011(H. 


H i - 1 . L Q.oiKH^KHt-DSTQ,; 
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where 
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1.268 
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0, e 
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— = 0. 123 
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-1. 561 H./M c n 0A'°- 268 /t\/- \0. 268 

e M e 0 M MV-M 

\ v 0 / \ t A m °/ 


(23) 


and x is the length along the surface of the body. 

Since the inviscid flow is established through the relaxative numerical calculations 
by repeatedly sweeping from upstream toward downstream in the flow field, it is expe- 
dient to account for the viscous flow effects by performing the boundary- layer calcula- 
tions at the end of each sweep (iteration) of the inviscid flow field with the prevailing in- 
viscid flow results on the equivalent body surface as the guiding free-stream information. 
A new body geometry is obtained through the displacement thickness correction proce- 
dure and is used for the next sweep of the inviscid flow field. The final flow pattern is 
established when the successive change of the disturbed potential is less than an arbi- 

c 

trarily small number (e. g. , 3. 0x10 ) for all grid points throughout the field. 
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METHOD OF CALCULATION 


The flow field in the g, £ -plane (0 < £ < 1, -2 < g < 3) is discretized into a system 
of finite points. There are 26 points in each column, which are uniformly spaced in the 
£ -direction. The suggestion of successive refinement of grid spacing by South and 
Jameson (ref. 6) is used for grid points along the £ -direction; thus, the number of 
grid points in each row changes from 26 to 51 and finally to 101. This distribution of 
points was used for all results shown in this report. The stretching parameter B 
(eq. (6)) was set to 8 for all results obtained from this series of calculations. Different 
values of B did not materially affect the results. Also, below a free-stream Mach 
number of 0. 9, there is no trouble in the convergence of the inviscid calculations toward 
the final solutions, and over-relaxation factors as high as 1. 5 were used to reduce the 
computational time. For the free-stream Mach number of 0. 9, some under-relaxation 
was required. However, because of the extremely sensitive character of the transonic 
flow, damping is required for the successive change in the equivalent inviscid body con- 
figuration when viscous flow equations are incorporated into calculations. This damping 
is more important at higher subsonic free-stream Mach numbers. A damping ratio of 
2:1 defines 


6 


* 

cor, i+1 


2 • 6 


cor 




(24) 


where 6* . i is the 6* to be applied for the correction to the body geometry for the 

forthcoming (i + 1) L inviscid calculation, 6* .• is the 6 added to the body geometry 

th * l 

for the i inviscid calculation, and 6* is the 6 calculated from the boundary -layer 
th ^ 

calculation on the i iteration. To illustrate the sensitivity of the results, the damp- 
ing and its effect, the intermediate variation of Cp values at the final point on the boat- 
tail (g = 1) is shown in figure 4. The results obtained with a damping ratio of 2:1 are 
shown in figure 4(a) for the fine mesh size. The pressure coefficient at the end of the 

xi_ 

boattail is diverging (rather than converging), and, on the 45 U1 iteration with the fine 
grid, program execution was terminated. When no damping was used in the calculations, 
termination occurred on the fourth iteration. No severe oscillations were observed in 
the coarse or medium mesh calculations, and convergence was relatively rapid whether 
or not damping was used. When the damping ratio was increased to 3:1 (fig. 4(b)), a 
slowly converging pattern was observed. Rapid convergence was obtained with a damp- 
ing ratio of 3:1 when boundary- layer calculations and displacement-thickness corrections 
to the actual body shape were made at every other iteration as shown in figure 4(a)). 

All viscous-flow calculations were based on the nominal flow conditions under which 
the model was tested. These conditions are given in the following table: 
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Free- 

Stagnation 

Stagnation 

Reynolds 

stream 

point pres- 

point tern- 

number 

Mach 

sure, 

perature, 

per meter, 

number, 

p o- 

N/m^ 

T 

0' 

Re/m 


K 


0. 56 

1. 24X10 5 

329 

1. 18xl0 7 

.70 

1. 19 



1.32 

.80 

1. 17 



1.39 

.90 

1. 15 



1. 45 


RESULTS OF CALCULATIONS 

The inviscid-flow results can be obtained by bypassing the viscous-flow calculations. 
The inviscid results for the four Mach numbers investigated are shown in figure 5, as 
are the results obtained when the viscous effects are taken into account. A substantial 
modification of the inviscid results due to the viscous-flow effects has occurred. In 
addition, the final results agree well with the experimental data, and the strong interac- 
tion character of the transonic flow past boattails is thus fully illustrated. For 

= 0- 9 a slightly higher pressure level in the last part of the boattail surface and an 
overshoot of shock recompression can be observed. 

The corresponding boundary-layer flow properties are presented in figures 6 to 8. 
Faster rates of increase of 6*, 9, H,, and H on the last part of the boattail surface 
are expected because of the adverse pressure gradient in this region. If ^ 2. 0 is 
taken to indicate flow separation, the results (fig. 8(a)) indicate that the flow has not 
separated from the boattail. 

To obtain these results (including viscous effects), the approximate amounts of com- 

— ft 

putational time required (with e = 3x10 ) for various free-stream Mach numbers M^ 

are 


Free- 

stream 

Mach 

number, 

Computer 
(CDC 6600 
CPU) time, 
min 

0. 56 

1.9 

.7 

2.0 

.8 

2. 1 

.9 

3.2 
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It has been observed that, when e is changed to lxlO - ®, the results are not signif- 
icantly different and it is not worthwhile to use computational time. On the other hand, 
when e is changed to 1x10“^, convergence may be indicated even when the oscillatory 
pattern has not been completely eliminated. 


CONCLUSIONS 

When viscous flow equations are incorporated into the calculations, damping the 
changes in the displacement thickness between iterations is definitely needed, particu- 
larly for cases with high subsonic free-stream Mach numbers. This is by no means in 
conflict with the conventional belief that the viscous flows are relatively insensitive to 
varying external flow conditions. It is the extremely sensitive character of the inviscid 
transonic flow that makes damping the change in the equivalent inviscid geometry 
necessary. In fact, for higher subsonic free-stream Mach numbers (M^ — 1), under- 
relaxation of the potential is needed (ref. 6) even for purely inviscid flow calculations 
(fixed body geometry). * 

At free-stream Mach 0. 8 and below, the theoretical boattail static -pressure distri- 
butions were in reasonable agreement with data, although the theoretical results indi- 
cated slightly higher pressures near the end of the boattail. At M,^ = 0. 9 this discrep- 
ancy is considerably larger. Whether the viscous flow analysis needs more refinement 
(such as normal pressure difference across the viscous layer or shock-wave - boundary- 
layer interaction) or the flow has been intermittently separated on this portion of the 
surface is yet to be examined. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, February 28, 1975, 

505-04. 
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APPENDIX - SYMBOLS 


AB1, AB2, AB3 
API, AP2, AP3 
B 

BBl, BB2, BB3 
BP1, BP2, BP3 


C 

C 

c ; 

C 

c f 

D 

D 


max 

D s 

E 

F 

H 

H c 

H i 

L 

M 

Pr 

P 

R 

Re/m 


function for boundary -layer calculations, defined in eqs. (23) 

coefficient functions at point where flow is subsonic, defined in 
eqs. (14) 

coefficient functions at point where flow is supersonic, defined in 
eqs. (15) 

radial stretch parameter 

coefficient functions at point where flow is subsonic, defined in 
eqs. (14) 

coefficient functions at point where flow is supersonic, defined in 
eqs. (15) 

dimensionless speed of sound 

1 2 

pressure coefficient, (p - p*,)/^ P^Voo 
pressure coefficient corresponding to local sonic flow 
local speed of sound 
skin friction coefficient 

coefficient in inviscid analysis, defined in eqs. (13) 
maximum model diameter 
sting diameter 

coefficient in inviscid analysis, defined in eqs. (13) 
coefficient in inviscid analysis, defined in eqs. (13) 
coefficient in inviscid analysis, defined in eqs. (13) 
shape factor of compressible boundary layer 
shape factor of incompressible boundary layer 
body length to end of boattail 
Mach number 
Prandtl number 
pressure 

boattail juncture radius of curvature 
Reynolds number per meter, P^V^/p,*, 
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r radial distance 

r^z) local body radius 
T temperature 

T reference temperature for the boundary -layer calculation 
U dimensionless velocity component in axial direction, u/V^ 
u velocity component in axial direction 

V dimensionless velocity component in radial direction, v/V M 
free stream velocity 

v velocity component in radial direction 

x distance along body surface 

z axial distance from upstream end of model 

y ratio of specific heats 

6* displacement thickness of the boundary layer 

6 * or correction to the body to account for viscous effects 
A? uniform spacing of the finite grid system in £ -direction 
A£ nonuniform spacing of the finite grid system in £ -direction 
£ transformed radial coordinate 

77 transformed coordinate, r - r^ 

0 , 9. momentum thickness of the compressible and the corresponding incompressible 
boundary layers 

p viscosity 

p viscosity evaluated at T temperature 

v kinematic viscosity 

| transformed axial coordinate 

<3? potential function 

cp normalized disturbance potential function 

Subscripts: 

e edge of boundary layer 

r partial derivative with respect to r 

w wall 
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z partial derivative with respect to z 

C partial derivative with respect to £ 

| partial derivative with respect to | 

0 stagnation state 
00 approaching free stream 
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joattailed afterbody 



Figure 1. - Geometry used for numerical calculations. 


O Data from ref. 1 


Small disturbance theory with boundary layer 



Axial distance, z/L 

Figure 2 . - Comparison of pressure distributions calculated 
from small disturbance theory with data at free-stream 
Mach 0. 80. 








Pressure coefficient, 




Axial distance, z/ L 


(c) Free-stream Mach number, 0.80. (d) Free-stream Mach number, 0.90. 

Figure 5. - Comparison of theoretical and experimental boattail pressure distributions. 
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Incompressible shape factor, 



Figure 6. - Calculated displacement thickness along the 
boattail. 



Figure 7. - Calculated momentum thickness along the 
boattail. 



Figure 8. - Theoretical shape factors along the boattail. 
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— National Aeronautics and Space Act of 1958 


NASA SCIENTIFIC AND TECHNICAL PUBLICATIONS 


TECHNICAL REPORTS: Scientific and 
technical information considered important, 
complete, and a lasting contribution to existing 
knowledge. 

TECHNICAL NOTES: Information less broad 
in scope but nevertheless of importance as a 
contribution to existing knowledge. 

TECHNICAL MEMORANDUMS: 

Information receiving limited distribution 
because of preliminary data, security classifica- 
tion, or other reasons. Also includes conference 
proceedings with either limited or unlimited 
distribution. 

CONTRACTOR REPORTS: Scientific and 
technical information generated under a NASA 
contract or grant and considered an important 
contribution to existing knowledge. 


TECHNICAL TRANSLATIONS : Information 
published in a foreign language considered 
to merit NASA distribution in English. 

SPECIAL PUBLICATIONS: Information 
derived from or of value to NASA activities. 
Publications include final reports of major 
projects, monographs, data compilations, 
handbooks, sourcebooks, and special 
bibliographies. 

TECHNOLOGY UTILIZATION 
PUBLICATIONS: Information on technology 
used by NASA that may be of particular 
interest in commercial and other non-aerospace 
applications. Publications include Tech Briefs, 
Technology Utilization Reports and 
Technology Surveys. 


Details on the availability of these publications may be obtained from: 

SCIENTIFIC AND TECHNICAL INFORMATION OFFICE 

ADMINISTRATION 


NATIONAL AERONAUTICS AND SPACE 

Washington, D.C. 20546 




